Oscillations in the Primordial Bispectrum I: Mode Expansion 
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We consider the presence of oscillations in the primordial bispectrum, inspired by three different 
cosmological models; features in the primordial potential, resonant type non-Gaussianities and de- 
viation from the standard Bunch Davies vacuum. In order to put constraints on their bispectra, a 
logical first step is to put these into factorized form which can be achieved via the recently proposed 
method of polynomial basis expansion on the tetrahedral domain. We investigate the viability of 
such an expansion for the oscillatory bispectra and find that one needs an increasing number of or- 
thonormal mode functions to achieve significant correlation between the expansion and the original 
spectrum as a function of their frequency. To reduce the number of modes required, we propose a 
basis consisting of Fourier functions orthonormalized on the tetrahedral domain. We show that the 
use of Fourier mode functions instead of polynomial mode functions can lead to the necessary factor- 
izability with the use of only 1 /5 of the total number of modes required to reconstruct the bispectra 
with polynomial mode functions. Moreover, from an observational perspective, the expansion has 
unique signatures depending on the orientation of the oscillation due to a resonance effect between 
the mode functions and the original spectrum. This effect opens the possibility to extract informa- 
tion about both the frequency of the bispectrum as well as its shape while considering only a limited 
number of modes. The resonance efi'ect is independent of the phase of the reconstructed bispectrum 
suggesting Fourier mode extraction could be an efficient way to detect oscillatory bispectra in the 
data. 



I. INTRODUCTION 



In recent years it has become evident that determining 
the precise physics of inflation requires the observation of 
higher order correlation functions beyond the power spec- 
trum [1]. These correlation functions can be obtained 
from the Cosmic Microwave Background (CMB) - 
12] and Large Scale Structure (LSS) |71 13, but recently 
[QHTT]. it has been shown that in principle 21-cm obser- 
vations of the early universe can also be used to measure 
n-point statistics. Because higher order correlation func- 
tions introduce more free parameters they can be used 
to constrain more complex models of inflation, since an 
increased set of parameters will allow for a unique fitting 
of the model to the observed data [T]. However, both 
due to computational and observational limitations, only 
the bispectrum has been reasonably investigated. For 
the detection of higher order correlations we will have to 
wait for more advanced data sets, such as Planck and im- 
proved analysis methods, although preliminary attempts 
have been made [3 [HI [13] . Even the detection of the bis- 
pectrum is not optimal, as a bispectrum would at least 
be a continuous three parameter observable but thus far 
only constraints have been set on limiting cases, in which 
2 of the parameters are fixed and the third one is mea- 
sured for a predetermined triangular configuration. The 
limiting cases (shapes) are known as the local, equilat- 
eral and orthogonal (and in the context of limiting tri- 
angular configurations; enfolded) non-Gaussian features. 
Precisely these features have been chosen, as it has been 
shown theoretically that most models of inflation produce 
non-Gaussianities that fall in one of these three classes 
(for recent reviews see P^HTO] ). 



When constraining non-Gaussianities using the bispec- 
trum, it has been a prerequisite that the comoving mo- 
mentum dependence should be factorizable; the bispec- 
trum should be separable into a product of functions of 
one variable, each variable being one of the three co- 
moving momenta making up the connected correlation 
triangle. Foremost, this requirement is set because of 
computational limitations that would render the analy- 
sis intractable if a given primordial bispectrum is not of 
the factorized form. The number integrals and sums one 
has to perform when computing an unfactorized bispec- 
trum scale with the number of pixels as iV^/^, while for 
factorizable shapes this reduces by one factor of N |17j . 
Although one integral can be computed fairly quickly 
the number of pixels (0(10^) for WMAP and 0(10'^) for 
Planck) is large and one factor of N can make all the dif- 
ference. The constrained bispectra, local, equilateral and 
orthogonal, have thus far been factorized templates. In 
case of equilateral [15] and orthogonal [Ij these have been 
constructed via approximation of a predicted signal, in 
the local case, the template is a direct representation of 
the theory [l9l - f2T| . For a particular type of bispectrum to 
be constrained, it is necessary to construct a factorized 
template that 'matches' the bispectrum. Until recent, 
there was no given prescription how to factorize a given 
theoretical bispectrum. In [22l [23] it was shown that 
factorizability can be achieved in both comoving momen- 
tum and multipole space by expanding the bispectrum in 
mode functions that are orthogonal on the domain of the 
bispectrum dictated by triangle constraints. The purpose 
of this factorization is to be able to quickly compute the 
full CMB bispectrum {Bi^i^i^) and generate CMB maps 
with a arbitrary primordial statistics (up to the trispec- 
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trum [13]) which are used to determine the variance of 
the (statistical) estimator. In the same paper, it was also 
shown that one can efficiently extract information about 
non-Gaussianity in the observed CMB by measuring the 
weight of each mode in the data and comparing this to 
theoretical predictions^. 

In this paper we investigate how well this mode ex- 
pansion works for a class of bispectra that contain (a 
large number of) oscillations. The reason to be interested 
in such features is that a number of theoretical models 
[H [IS [371 US ED predict oscillations in the bispectrum 
and in order to be able to constrain such models, a plau- 
sible first step is to factorize these bispectra. As it is, os- 
cillations can be considered as an extra, distinguishable, 
degree of freedom within the bispectrum which could re- 
sult in narrowing down the number of potential scenarios 
of infiation. 

We introduce three different cosmological scenarios in 
which oscillations in the bispectrum can appear. We will 
briefly discuss the theory behind these models and show 
to what extend these would be distinguishable from one 
another in the data in section |lll Two out of three bis- 
pectra can have significant correlation and it could be 
difficult to discriminate between such models in future 
surveys. We will discuss the method of polynomial ex- 
pansion in order to rewrite the primordial bispectra in 



factorized/separable form in section III As expected, the 
number of modes required in the expansion grows along 
with the frequency of the theoretical spectra. In sec- 



tion III A we show how fast polynomial expansion would 
yield a reasonable reconstruction of the given bispectra 
predicted by the three cosmological scenarios. Subse- 
quently we will investigate another set of modes that can 
lead to a separable expansion of the theoretical bispec- 
trum in section IIIIBI These modes are based on the sine 
and cosine and the resulting set of orthonormal functions 
can be considered a Fourier-type basis on the tetrahedral 
domain. After detailing the construction of this set of 
orthonormal mode functions, we will compare the num- 
ber of modes required to achieve comparable correlation 
with the polynomial mode expansion . It turns out that 
this number is reduced significantly and as such Fourier 
expansion can be considered a reasonable alternative to 
expand oscillatory spectra. For larger frequencies both 
Fourier and Polynomial mode expansion become ineffi- 
cient. Fortunately, it turns out that for various oscilla- 
tory signals only a limited number of modes contribute 
significantly in the reconstruction of the original spec- 
trum. This has various consequences for the viability of 
Fourier mode expansion as well as possible observational 
advantages compared to polynomial modes, which will 



^ During the finalization of this paper, the same group published a 
paper I31| in which many non-factorizable non-Gaussian shapes 
have been constrained using the WMAP 5 year data and the 
method of mode expansion 



be discussed in section IV We conclude this paper in 
section IVl 



II. OSCILLATIONS IN PRIMORDIAL 
BISPECTRA 



In this section we will briefly discuss 3 distinct possi- 
bilities that can produce non-Gaussianities that have an 
oscillatory component. Two of these examples have an 
exact solution, while a third has only been solved numeri- 
cally and we will use an approximate form. In the follow- 
ing paragraphs we will describe the physics behind these 
models and quote their theoretically predicted primor- 
dial bispectra. In addition we investigate how well these 
bispectra can be distinguished from one another by com- 
puting their correlation, which will be defined shortly. 
Since all these bispectra have poor overlap with existing 
spectra, there exists substantial room for improvement, 
which we could achieve by approximating these shapes 
via mode expansion. This will be the topic of the next 
section. 

For completeness, let us introduce (standard) notation. 
The primordial bispectrum is given by 

(1) 

where ^ is the gauge invariant curvature perturbation 
(C = —H6(j)/(j)o) which is constant after horizon exit, A 
is the amplitude of the primordial power spectrum (i.e. 
for single field slow-roll A = H'^/Sne, where H is the 
Hubble rate at the end of inflation and e the slow-roll 
parameter) and F(ki,k2,k^) is the shape of the bispec- 
trum. We will also make use of 5 = k\k\k'^F. In the 
following we will discuss the shapes of the bispectra and 
quote theoretically predicted ranges of their associated 
Jnl- We would like to refer to the literature for a de- 
tailed examination of the theoretically predicted values 
of /at l ^24, ,25, i271 HSj in various theoretical contexts. 



A. Features in the Potential 

Sharp features in the potential can temporarily break 
slow-roll and produce large non-Gaussianities [26l[27]. As 
long as the system relaxes within several Hubble times, 
inflation can still lead to a significant amount of e-folds to 
solve the standard cosmological problems. The motiva- 
tion for these type of features is two-fold. First, there are 
hints of glitches in the primordial power spectrum that 
could be cross-checked using the bispectrum [3S]. A sec- 
ond motivation is theoretical in nature. In certain brane 
inflation models the effective 4-dimensional potential dis- 
plays sharp features (see [H] and references therein). 
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One of the possible sharp features is a step in the po- 
tential, which can be parameterized as 



v{4>) 



1 
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1 



; tanh 



(2) 



where c, d and (ps respectively determines the height, 
width and location of the feature. 

The resulting bispectrum can only be computed nu- 
merically. The authors of ^27) have proposed an approx- 
imate analytic form 
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(3) 



The approximation can in principle be improved |14| 
by multiplying by an 'envelope' function, but such im- 
provement would not gain us any more useful insight re- 
quired for the analysis in this paper and we will there- 
fore omit it. Here /c* is related the location of the fea- 
ture in the potential (j)s- Evidence for features in the 
power spectrum around I ^ 30 have been put forward 
in [3S]. It was shown that the inclusion of features 
in the primordial potential could improve the best- 
fit. Such a feature would approximately correspond to 
fc* = 30/?7o ^ 0.002Mpc~"'^. This relation also indicates 
that the smaller the scale at which the feature appears 
the larger the associated wavelength. Roughly the wave- 
length corresponds to the location of the feature, e.g. for 
a feature at / = 30 the wavelength SI ~ 30 [3H]. Here 
we do not necessarily relate to an observed feature at a 
specific value in multipole space since features that lead 
to non- vanishing bispectra can still be present with min- 
imal consequences for the observable power spectrum. 
The quantities we will compute in the remainder of this 
paper are mostly integrals that run over the domain of 
comoving momentum space between kmin k < kmax- 
It is therefore convenient to choose our reference scale 
kmax ~ 10~^Mpc^^, the smallest observable scale in the 
data, in order to be able to compare the frequencies 
in the various models. We then define xi — ki/kmax, 
X2 = k2/kmax, X3 = k^/kmax, Xt = h / kmax and rewrite 
the shape of this bispectrum as 
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k: 



_g sin{ujfXt + S) 



999 



(4) 



with oj-f — kmax/k*. For a feature at k^, — 0.002Mpc^^ 
we therefore find w/ ^ 50. Note that lu — 50 can be 
considered an upper limit in allowable frequencies due to 
features in the potential. For features at smaller scales 
the frequency will be smaller. This bispectrum with a 
frequency of w/ = 50 is shown in the bottom of figure [T] 
The amplitude of this type of non-Gaussianity is gov- 
erned by the width and the depth of the feature in the 
potential 
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de ' 



(5) 



which for a feature at Z ^ 30 would imply /j^*^' ~ 0(10) 



B. Resonant non-Gaussianity 

This type of non-Gaussianity is a result of a periodic 
feature in the inflaton potential as apposed to a sharp 
feature explored in the previous example. These features 
will cause oscillations in the coupling(s) of the interaction 
terms of the inflaton field. Resonance occurs when an 
oscillatory mode well within the horizon grows during 
inflation until its frequency hits the same frequency as 
those of the couplings. So as long asw > H resonance will 
occur at some point within the inflationary history of the 
mode. This resonance can result in a large contribution 
to the three point correlation function [57] . 

In a general scenario, with an oscillatory potential we 
obtain an expression for the bispectrum of the form |271 



J 7-, 



4-2 i.2 1,2 



sin(Cln(fcf/fc*)) 



+C-^cos{Cln{kt/k,))J2'^ 



(6) 



Here C is related to the frequency a.s C — uj/H with H 
the Hubble rate during inflation (which is approximately 
constant) and fc* introduces a phase. One can also com- 
pute the general expected amplitude of non-Gaussianity 
which is related to the frequency as 
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(7) 



Here rjA represents the amplitude of the oscillatory com- 
ponent of the couplings. 

Physically such features might be realized in terms of 
brane inflation [35] where the periodic feature comes from 
a duality cascade in the warped throat, as well as axion- 
monodromy inflation where the periodic feature is a re- 
sult of instanton effects [28l [30] . As an example let us 
consider the latter. Axion inflation is well embedded in 
string theory and represents a favorable candidate for in- 
flation if the observed tensor modes are relatively large 
(r ~ 0.07). Such a scenario implies inflation occurred at 
energies close to the GUT scale and would indicate that 
we require the knowledge of the UV completion. 

The axion potential is given by 



y(0) = T/o(0)+A'cos0//. 



(8) 



The parameter / represents the axion decay parameter. 
The range of / which would generate observable non- 
Gaussianities and is still consistent with observations of 
the power spectrum is given by 10"'* ^ / ^ 6 x 10~^ 
[IHl [30] • The lower bound is set by the requirement that 
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the period of the oscillation should be larger than Al 1 
for I = 200. For a linear zero order potential the resulting 
bispectrum is then given by 



sin(a;r Inxt + 7i) 



+UJr ^ C08{uJr InXf + 7l) ^ — 



(9) 



with (jjr — {f'f>*)~ and 71 = lj,^. In kmax/k*. Here k^. is 
pivot scale (/c* = 0.002Mpc^^), and 0* is the value of 
the inflaton field when the pivot scale exits the horizon 
and is of order 10 (Mp). Given the range for /0* the 
frequency of the oscillations in the bispectrum lie within 
20 < lO'^. A plot of this shape is shown in the top 
right figure [1] 

The amplitude of the axion bispectrum (for a linear 
potential) is given by 
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(10) 



The amplitude is therefore proportional to a power of the 
frequency. For a linear potential b — A^/(/i'^/), where 
H ^ 6 X 10~^ is fixed by COBE normalization. From 
observations of the power spectrum one can constrain 
bf < 10-^ 'W^ and therefore 



fres 1 n-3 5/2 

JnL ^ 



(11) 



allowing 0(1) < f^'^l < O{l0^). 



C. Initial State Modifications 

Since inflation is an effective field theory in a curved 
background, choosing an appropriate vacuum state is by 
no means evident^. In general the initial or vacuum 
state is chosen to be equivalent to the free field vac- 
uum state in flat Minkowski space, know as the Bunch 
Davies (BD) vacuum. Although it seems that possi- 
ble corrections to this assumption are constrained to be 
small (from general observation of the power spectrum 
[551 [5^ and backreaction constraints [Mllll])i it has 
been shown that small corrections in the BD state can 
result in rather large non-Gaussian effects [2^ E51 [Ml 157] . 
Using the currently available bounds on non-Gaussianity 
from CMB data, deviations from a pure Bunch Davies 
stae have been constraint even further, although these 
constraints strongly depend on the inflationary model. 
However, there exist significant room for improvement 



^ For a in depth discussion on deviations from a BD vacuum see 
for example |25l I41| and |14| section 6.4 



as non-Gaussianities from these modifications are highly 
oscillatory and therefore the derived constraints are rel- 
atively poor since they depend on the correlation with 
measured smooth bispectra. 

A number of different scenarios have been considered in 
which initial state modifications were investigated. Here 
we will not discuss all of these, although the results can 
differ significantly Such differences make it difficult 
to make robust predictions, it seems inevitable however 
that once you introduce a effective field theory cutoff, 
oscillations appear in both the power and bispectruni. 
We will consider one example that represents a large 
class of models with a non-canonical effective field the- 
ory action, which already drives large non-Gaussianities 
to start with. This particular class has a speed of sound 
Cs < 1, such that perturbations in the medium propagate 
slow compared to the growth of the causal horizon. The 
leading order shape of the resulting bispectrum is given 
by [25j 



F, 
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1 cos{kjCsr/o + S) 



fclfc2^3 \ 2 

sin(fcjCs77o -t- S) 
(kjCsrioY 



kjCsTJo 

cos S — cos(fcjCs77o + 6) 
[kjCs-qoY 



(12) 



Here kj — kt — 2ki. In 124, it was assumed that there ex- 



ists a fixed physical cutoff hyper-surface rjo that is scale 
dependent such that the overall momentum dependence 
of the bispectrum becomes scale invariant. Such a choice 
is known as the New Physics Hypersurface (NPH), as ap- 
posed to Boundary Effective Field Theory (BEFT) ap- 
proach in which the cutoff is time dependent [40]. The 
subtlety is that the cutoff appears due to the presence 
of a non-BD state in each direction in comoving momen- 
tum space. Consequently, ??o(fci) will depend on the ki 
direction the BD vacuum has been perturbed in. This 
direction is set by the direction in which ki picks up a 
minus sign due to the Bunch Davies vacuum perturbation 
as explained in |24' . One could allow for scale invariance 
breaking and consider BEFT, however there are some 
suggestions [35] that such large scale invariance should 
have been observed already. We can rewrite the bispec- 
trum as 
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^2 (^.±1±^^ _ 1) 

i<5-C0S (,.„^2±l±^^+72) 

u;3(^i±i±^-l)' 



(13) 



where 72 = 6 — tUy and ujy = krjQCs — {k / a^) / {H / Cg) or 
the ratio between the largest physical momentum scale 




FIG. 1: 3 examples of oscillating bispectra. We have set u)^ = ujf = cjr = 50. The pivot scale in eq. ([9| is set to fc» = 0.002Mpc~ 
and a;i = 1 — >■ fci = kmax- In addition 5 = in both the modified initial state and feature scenario. The non-BD bispectrum 
contains the most features and, not clear from this image, the number of features (effective frequency) increases rapidly for 
smaller values of xi making this bispectrum particularly hard to reconstruct using mode expansion. 



and the Hubble radius at time rjQ which can be as large 
as 10^ [231 113 131] ■ Note that from this expression it seem 
that Xj = Xj+i +Xj+2 represents a singular line (the en- 
folded limit). However, one can show that all infinities 
are cancelled against each other and the the expression is 
finite and vansihing'^.When computing quantities numer- 



^ This limit is on the enfolded line Xj = Xj^i + Xj^2 within the 
sum. Outside the sum, this expression is non-zero but finite. For 
example xi X2 + x^ gives: 

(("^ + - ' (("^ + - cos (^^) 

/, 2 2 2\ f "ix^iUJ \ 

-X2 [{X2 + x-i) - 2x-j^u] ) COS 

V \X2+X-iJ 

, / / 2X2U \ ( 2X3UJ WW 
+ 2x3{x2 + Xijuj sm + sm ■ 

V \X2+X3J \X2+X-iJJJJ 



ically, such as the correlator in section 3, these apparent 
singularities can be hard to handle and we need to be 
aware of these. We have plotted this shape in the top 
left figure [T] 

The amplitude of the non-BD bispectrum is a function 
of the frequency and the Bogoliubov parameter quanti- 
fying the deformation away from the BD state. The way 
this bispectrum was computed, considered a Bogolyubov 
correction of linear order /3 and small speed of sound . 
In this particular scenario, /nl is roughly given by 



f^l^ ^ ^ulp. (14) 



From backreaction and power spectrum constraints /? ;< 
10~^, which could still allow observable levels of non- 
Gaussianity. 
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D. Distinguishability 

Although the presented theoretical bispectra have dif- 
ferent characteristics, we would like to get an indication 
how well these could be discriminated. For instance, it 
seems obvious that the similarity between the feature 
bispectrum and the resonant bispectrum could lead to 
significant confusion when actually traced in the data. 
In order to do so, we want to measure the distinguisha- 
bility of these shapes, which is usually quantified using 
the amount of overlap or correlation between two shapes. 
One can define a inner product between two shapes 



Fx * Fy = I dkidk2dk^k\k2k^WkFxFY 
'A. 



dkidk2dk^WkSxSy 



(15) 



The correlation between two shapes Fx and Fy is then 
defined as 



G{Fx,Fy) 



Fx * Fy 



(FxirFxy^^{Fy*Fy)y^- 



(16) 



Here Wk is a weight function, which was chosen as Wk = 
1/kt in [23j to increase resemblance with the Fisher ma- 
trix (correlation) found in multipole space. The integral 
runs over the 'tetrahedral' domain, which is bounded by 
the following triangle constraints 

ka <h + kc for ka > fcb, kc 

kaj ^6; k^ ^ kjyif^xt 

where a,b,c = {1, 2, 3}, a ^ b ^ c. 

Before we compute the correlation between the shapes, 
let us perform a quick qualitative analysis in order to get 
an indication of what to expect. First of all, note that the 
shape coming from initial state modifications (cq. (13)) 
is clearly different from the other two. While for features 
(eq. Q and eq. ([9])) the argument in the oscillating 
functions explicitly depends on the sum all three comov- 
ing momenta, the argument in eq. ( 13 ) depends on the 
ratio of momenta. Consequently we can expect a rather 
small overlap. This becomes even more apparent once 
we adapt a new set of variables 

k^kt/2, ki^k{l~P) 



k2 = -fc(l + a + /3) 
dkidk2dk^ 



kz = h{l-a + l3) 
k'^dkdadp, 



proposed in |22j . As a consequence the argument in eq. 
( 13 1 will depend on the two variables a and /3, while the 
arguments in eq. Q and (|9| will only depend on k. In 
that sense, we can say that oscillations in these shapes 
are in orthogonal directions. 

In addition, for both the feature and resonant bispec- 
trum the frequency is fixed along one direction. That is. 



the frequency does not change (feature) or only slightly 
changes (resonant) when you run through a fixed direc- 
tion in comoving momentum space. For the non-BD bis- 
pectrum however the argument in the oscillating function 
has a component that scales as l/xj. Consequently for 
the effective frequency uJeff — ^ oo- Naturally, 



is cutoff from below (as kmin/k„ 



10 ), however 



even with a cutoff the range in effective frequencies is 
large along a direction. This effect is present at all fre- 
quencies, and it turns out it will determine the efficiency 
of mode expansion for this bispectrum discussed in the 
next section. 

We have numerically calculated the correlator as de- 
fined in eq. ( 16 1 between both feature bispectra and 



non-BD spectrum. We found the correlation to be maxi- 
mal for low values of both frequencies (of order 1 percent 
around uj — 10), indicating that there is no evidence for 
a particular resonant frequency; the largest correlation 
occurs due to the fact that there are less oscillations, 
thereby decreasing the chance for (almost perfect) can- 
celations in the integral. As expected, we can safely con- 
clude that these shapes are distinguishable/orthogonal. 

For the two bispectra of eq. Q and (|9| we can expect 
a larger correlation. The appearance of a log in eq. ([9]) 
is the only major difference between the two bispectra. 
In the new coordinate set, the bispectrum of Q does 
not depend on the a or /3. Let us try to make a simple 
analytical approximation of the relevant correlator before 
we compute the correlation numerically. The first term 
in ([9| dominates the second for large values of (/0*)~^. 
Therefore for simplicity we neglect the second term. As 
a consequence both terms now depend only on k. In the 
computation of the correlator the integration over a and 
/3 drops out and to get an indication of the resonance we 
only need to investigate the following integral: 



Xtdxt sm{ujfXt + S) sin(a;r log xt + 71). (17) 



where we assumed that at most kf — ^^k^^x — ^ Xf — 3. 
This integral can be done analytically and results in a 
sum of r functions (we have set (5 = 71 =0). The inter- 
pretation of the result is rather complicated as all terms 
are divergent and there are no terms that can be easily 
neglected. However, one can plot the result and find that 
there is a clear resonance 'area' around lUj. ~ 20^^. We 
have confirmed this resonance as a function of frequency 
when considering the full expression and allowing both 
phases to be non-zero. We have plotted (fig. [2]) the cor- 
relation for a range of frequencies {10 < u < 1000) and 
a phase (5 = 71 =0. The largest values obtained from 
this numerical computation are of order 0.6, or 60 per- 
cent correlation (we have used discreet steps of Sut = 10), 
and we expect there to exist correlation of 0(1) for some 
specific values of w). As such it will be hard to discrimi- 
nate between these two models solely using observations 
of the bispectrum (as one could simply confuse frequen- 
cies). However, as mentioned before, axion inflation for 
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FIG. 2: The correlation between the bispectral shapes of Q 
and ([9| for various values of the frequency. The light ar- 
eas correspond to correlations of order 0{1), while the dark 
shaded areas correspond to correlations close to 0. The cor- 
relation was computed with 5 = 71 = 0. 



example predict a large scalar to tensor ratio. Measure- 
ment of r could break the degeneracy between a sharp 
feature in the potential versus axion inflation. In addi- 
tion, one does not expect uJr < 10 since it would not 
produce observational f]^l, while for the feature bispec- 
trum the natural frequency is no larger than ujf ^ 50. If 
one would be able to extract a frequency from the data, 
a large frequency would favor a resonant model while a 
low frequency could indicate a sharp feature. 



III. MODE EXPANSION 



A. Power Modes 



indicated bispectra with completely orthogonal charac- 
teristics. Of course, optimally, one would simply look for 
the full bispectrum as a function of the multipole num- 
bers instead of constraining the amplitude in particular 
bispectral configuration, but the low S/N and computa- 
tional limitations have so-far restrained us to the former. 

The second reason not to look for more 'exotic' bis- 
pectra is that for a fast estimator, the bispectrum one 
would like to constrain needs to be factorizable and scale 
invariant. That is, it is useful if the bispectrum can be 
we written as sum of products of functions, where each 
function only depends on one direction in multipole or co- 
moving momentum space. It has been shown that such 
factorizability reduced the number of computations one 
has to make in order to constrain the amplitude of the 
bispectrum by a factor P, where / is the number of ob- 
servable multipoles of the experiment (leaving only P 
computations). 

The constrained non-Gaussian amplitudes (in the form 
of /^i, where i labels the comoving momentum type, lo- 
cal, equilateral or orthogonal) are all based on templates 
that are factorized in the manner explained above. For 
instance, although DBI inflation docs not produce a fac- 
torized bispectrum, it is well approximated by the equi- 
lateral template [18] , that is factorized by construction. 
The same is true for both the local and orthogonal tem- 
plate, as well as the enfolded [23] template. However, the 
method for constructing such factorized approximations 
of existing theoretical bispectra is rather ad- hoc. Until 
recently there was no procedure no construct a factorized 
bispectrum using a consistent prescription. 

In |23| . a method for constructing factorized approxi- 
mations to theoretical bispectra has been proposed using 
polynomial expansion. The approach is fairly straight- 
forward; one defines a set of orthonormal 3 dimensional 
functions (where orthonormal is defined using a correla- 



tor of the form eq. (16), and the weight function can 



be adjusted) which are a-priori factorized and from there 
one computes the corresponding weight factors (a„) via 
the inner product between a number of polynomial modes 
(Rn) up until a sufficient overlap between the polynomial 
expansion and the original bispectrum is established, i.e. 
until N such that 



The discussed primordial bispectra have very little in 
common with the constrained local, equilateral and or- 
thogonal bispectra. Typically, to constrain any type of 
non-Gaussianity one computes the correlator (eq. (16)) 



and derive the so-called 'fudge' factor which indication 
how much 'signal' leaks into an existing template With 
the use of the fudge factor one is able to deduce a bound 
on the amplitude of the unconstrained bispectrum. The 
reason why certain templates have been constrained and 
some others have not, is two-fold. First and foremost, 
until now most models produced non-Gaussianities that 
can roughly be placed in one of the constrained types. 
For this reason, it was not immediate to search for any 
other type, simply because there were no models that 



S{xi,X2,X3) 



N 

E 

n=0 



a„i?„(xi,X2,X3). 



(18) 



* For the construction of these polynomial modes we set w = I. 
Once computing the correlator between the original and the re- 
constructed spectrum one can take w = 1/kt in order to see how 
much of an effect projection onto multipole space can have. We 
find that it reduces the correlation by 5 to 10% in both polyno- 
mial expansion and Fourier expansion. As such, it should not 
effect the conclusions we draw in this paper where all correlation 
shown are based on w = 1. In order to build modes that are op- 
timized for multipole expansion you should start by considering 
a weight function 1/fet. This is beyond the scope of this paper. 
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Without discussing the details of constructing such poly- 
nomial modes (see for a detailed description), here 
we want to try and investigate how well this would work 
in case of oscillatory bispectra of Q, ([9| and (13 1. 

Before we do so, let us make a few notes. First of all, 
recall that the objective of the expansion is to factorize 
a given theoretical bispectrum. However, as you can see 
from eq. Q , this particular bispectrum, albeit a best-fit 
approximation^, is already of the factorized form. One 
can still try to expand this in terms of power law polyno- 
mials, as described here, since polynomial modes will in 
general behave better numerically. The other two exam- 
ples of primordial bispectra are not factorizable in terms 
of oscillating functions using simple identities. Conse- 
quently, the polynomial expansion seems to be a good 
first effort in order to set up an approximately factorized 
form. 

Secondly, were we able to expand these into a factor- 
ized form, and subsequently projected to multipole space 
and applied to the data, we might still miss the entire 
signal, simply because one of the free parameters is the 
frequency of the oscillations. For a non-BD bispectrum 
and the axion inflation model, the range of possible fre- 
quencies spans (at least) 2 orders of magnitude. There- 
fore, if we would fix the frequency, searching for a signal 
with a constructed factorized template would probably 
not be the best approach. Fortunately, we will later see 
that if you would measure mode functions in the data, 
instead of a fixed template, one could in principle extract 
information about a variety of oscillating signals. Let us 
emphasize that even if we would not be able to recon- 
struct a factorized form of a given spectrum with a small 
number of modes, it is still very well possible we could 
observe the same spectra by measuring a small number 
of mode functions in the data (effectively the frequency 
(and the phase) remain a free parameter during mode 
extraction). 
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2 


3 


4 


5 


6 


7 


8 


9* 


# of modes _R„ 


1 


5 


8 


12 


18 


43 


55 


69 


82 



TABLE I: As the frequency is increased it requires a rapidly 
growing number of modes to get over 98% correlation with 
the original spectrum. 



table |T] we have computed the number of modes neces- 
sary to get a correlation of at least 98% with the original 
spectrum for several values of w/. As expected, as the 
frequency is increased, one has to expand the bispectrum 
with a (rapidly) growing number of modes. For ujf — 9 
we get a 93% correlation with 82 modes. On itself, it 
actually quite remarkable that one is able to reproduce 
the spectrum with a limited number of modes. Recall 
that the possible feature at Z ~ 30 would result in a (de- 
caying) oscillation with ujf ^ 50, which would be hard to 
fit this way. On the other hand, as we argued earlier, a 
frequency oi iOf — 50 can be considered an upper limit, 
as features at higher multipole number would result in 
longer wavelengths. 



2. Resonant Bispectrum 

Next, let us consider the resonant bispectrum. It is 
quite similar to the feature bispectrum, but theoretically 
we expect much larger frequencies (20 < ujr < 10'^). 
We have computed (figure |4]) the correlation between 
the expansion and the original spectrum, chosen to be 
sm{uJr Ina;*) since for the same reason as before the phase 
will barely affect the number of modes required to recon- 
struct the spectrum. As expected, the convergence of the 
correlation towards one (perfect overlap) proceeds slowly. 



1. Feature Bispectrum 

First we consider the bispectrum coming from a fea- 
ture in the potential (eq. Q). Out of the given exam- 
ples it has the simplest form (excluding the envelope). 
We choose 5 — for simplicity, and since the phase can 
always be scaled out it will not affect the results^. In 



The proposed envelop function has the form (fci +^2-1- 
I~^^n^(ki+k2+ks)/ktm where m and n are fitted to the numeri- 
cal results. The envelope function is therefore also factorizable. 
Again, we did not consider this envelop since it is smooth com- 
pared to the oscillatory part of the bispectrum. However, such 
an envelope could be of significant influence in predicting the 
correlation in multipole space [22]. 

One would also have to consider cosujfXt but wc found no dif- 
ference when expanding between the cosine and sine in terms of 
the required number of modes. 




20 40 60 

# of modes Rn 



FIG. 3: Example of the increasing correlation (eq. ( 16 1 be- 
tween the approximation of eq. Q and the original spectrum. 
Here uuf = 9 and we find that it requires over 80 modes to 
achieve perfect correlation. 
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FIG. 4: The correlation between sin(a;r Init) for 3 diflFerent 
frequencies. From top to bottom uir = 20, 40 and 60. Be- 
yond frequencies of 60 polynomial expansion would require 
many modes to achieve significant correlation with the origi- 
nal spectrum. 
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FIG. 5: The correlation between the non-BD bispectrum of 



eq. ( 13 \ and its polynomial expansion as a function of mode 



number for 3 different frequencies aj„ = 20 (solid), 40 (dashed) 
and 60 (dot-dashed) . We have set 72 = but have found very 
little difference for non-zero 72. 



For the lowest frequency we considered (w^ — 20), the 
correlation reaches 71% after 82 modes. For LUr — 60 the 
largest correlation we can achieve is 7% after 82 modes. 
Recall that the amplitude of the resonant bispectrum is 
proportional to its frequency. The maximum correlation 
between existing templates and the axion spectrum is of 
order 1% ^30 (although for small frequencies this can be 
10% for the equilateral template) and possibly measur- 
ing these modes in the data would therefore still allow 



for a constraint on axion inflation that is 10-100 times^ 
better than what we have now. In general, increasing the 
frequency above uJr ~ 60, large correlation becomes hard 
to achieve with a limited number of modes. 



3. Non-BD Bispectrum 



For the non-BD spectrum of eq. ( 13 1 we find that com- 



puting the correlator numerically requires a very high 
resolution, because this shape contains terms that are 
singular and the spectrum as a whole is only finite due 
to the exact cancellation between all the specific terms. 
To avoid these problems one has to stay away from the 
line(s) Xj+i -I- Xj+2 = Xj, which can be done by adding 
a small e in the vicinity of this line in the integral that 
defines the dot product (eq. ( 15 1). The results are shown 
in figure [5] Even for low frequency {ojy = 20) we can not 
achieve a large correlation with 82 modes. On the other 
hand, increasing the frequency does not really affect the 
ability to reach similar correlation. Overall, we find that 
the non-BD bispectrum is the most difficult to recon- 
struct due to the appearance of terms that diverge inside 
the argument, since for Xj — > the frequency of the signal 
becomes extremely large at some of the edges of the tetra- 
hedral domain. The observation that we can still reach 
some correlation is because there are also areas on the 
tetrahedral domain where the effective frequency is rela- 
tively small. These areas remain even if ujy increases (al 
though they should become smaller and smaller) Conse- 
quently we find that achievable correlation with 82 modes 
is small but does not decrease significantly when you go 
to higher frequencies. The correlation with smooth spec- 
tra is typically of order ^1% [25j which means that an 
accumulated correlation of only a few percent could dras- 
tically improve the constraints we can put on as the 
amplitude scales proportionally to uj^. 

Polynomial expansion seems to work reasonably well 
for low frequencies of the various bispectra. For larger 
frequencies, to reconstruct the original spectrum the 
polynomial expansion requires an increasing number of 
modes. Given the large allowable frequencies for the reso- 
nant and non-BD bispectra, polynomial expansion might 
not be the most effective way of expanding. In the next 
section we will explore another type of expansion which 
uses a Fourier basis. We will investigate if such a basis 
would require less modes to achieve similar correlation. 



^ This would require adding the modes once extracted from the 
data. 
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TABLE II: The association of mode numbers for the first 27 modes. For example n = 10 corresponds to the mode for which 
one direction is of maximally 4th order and the other two are order, i.e. for polynomial modes -Rio oc f{l,xi,Xi,Xi,Xi) + 

/(I, X2,xl,X2,xi) + /(I, X3,X3,X3, xj) . 



B. Fourier Modes 



The polynomial expansion of j23j is based on power 
modes, i.e. the expansion is in increasing order of a;". 
This is not necessarily optimal for describing oscillatory 
functions. There are two possible alternatives; the first 
one would be to expand the argument into a sum of func- 
tions, that each depend on one direction only. A such, 
one can again use trigonometric identities to expand the 
cosine and sine into factorized forms (be that oscillatory 
functions). The second option could be to use a Fourier 
expansion instead of a polynomial expansion. This would 
only be useful if for large frequencies you would need a 
small(er) number of modes. Before we get into Fourier 
mode expansion let us briefly discuss the alternative of 
expanding the argument in the oscillatory function. 

This option would only suffice if the approximation re- 
quires 2 modes maximally. If it requires more modes, 
you will get product of two or three different directions 
in momentum space, and as a result you will not be able 
to expand the cosine and the sine. Let us consider the ax- 
ion model. The argument is given by ujlogxt- Using the 
mode expansion, one finds that one can achieve > 99% 
correlation after just two polynomial modes; zero order 
and first order. Not surprisingly this is almost equiv- 
alent to a Taylor expansion to first order of logfcj/fc* 
around the point k-t ^ 1.4 ~ ^/2. Consequently, there 
are no cross-terms, and one can expand the cosine and 
sine into factorizable function of the three comoving mo- 
menta, just like you could expand the feature spectrum 
into oscillating functions. As it turns out however, al- 
though there is a 99% correlation between the arguments 
after expansion, the full bispectrum is very sensitive to 
small deviations in the argument, especially for large fre- 
quency. Consequently, the correlation between the full 
bispectrum and the approximated bispectrum decreases 
as a function of the frequency; from ^ 90% for LOr = 1 
to ~ 50% for u!r = 20. Although this is equivalent to 
what can be achieved with the polynomial expansion us- 
ing 7 modes, the problem is that we can not improve 
it in any way. Since this will only work for a first order 
expansion, we can never reach beyond 50% correlation, 
unlike the polynomial expansion, where we can simply 
include more modes. Note that for non-BD model this 
method will not work as the argument is already a prod- 
uct of two directions in comoving momentum space, i.e. 
{kj+i + kj+2)/kj. 



The second option is to consider a Fourier expansion, 
where we try and fit terms such as sin ujf{x, y, z) to a sum 
of Fourier modes that all depend on one direction only. 
Such factorization would still lead to the P reduction 
in computation, since the integrals in k space can now 
be performed individually*. We consider exp[i27mx] as 
our basis function (as apposed to x") and constructed a 
orthogonal set of three dimensional mode function similar 
to [231 . The first few one dimensional functions are given 
by 



fo{xi) 

f2{xi) 

h{xi) 



= 72 



0.22 + Q.Tii + lAbe^"''^ 

-0.0087 + 0.041i(0.088 + 0.62i)e^'''"^i' - 

(0.31 + l.l2i)e^"''^ 

(-4.9 4- 2.2i)10~^ - (0.15 - 0.lli)e'^^''''^ -|- 
(0.68 - 0.65i)e''*''^i - (0.59 - 0.65i)e^"="i 
(-6.5 - 3.8i)10"'' - (0.042 + Q.{)m)e'^"''^ 



(0.44-f 0.12i)e^ 



(l.-H0.2i)e^"''i + 



(0.63-f 0.1i)e*^'^^i 

(-1.1 - ll.li)10~^ - (0.002 + 0.011i)e2"='i) 
(0.051 -I- 0.16i)e''"^i - (0.25 + QMi)e^''"'^ + 
(0.4 + 0.89i)e*^'^"=i - (0.21 -I- 0.4H)ei°"^i 



The functions are shown up to n = 10 in figure [6] From 
these one can construct the three dimensional basis func- 
tions via a product of each mode and symmetrization of 
three comoving momentum arguments; xi, X2 and 

Zprs{xi, X2,X3) cx [/p(xi)/r(x2)/s(x3) + 5 perm]. 

One has to introduce a counting scheme to re-numerate 
the three labels {p, r, s} to n. We have chosen equal slic- 
ing counting [3T|, of which the first 27 modes (n) and 
their association ({p, r, s}) are shown in table |ll] After 
the construction of these modes, one has to apply ad- 
ditional Gramm Schmidt orthogonalization to to in- 
crease orthonormality of different mode functions. We 



* In '23' Fourier mode expansion is briefly discussed in section E 
as a possible orthonormal basis, liowever no results are shown. 
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can be rewritten into a product of Fourier modes simply 
by using trigonometric identities, e.g. 

sina;/a;t = coscjyxa (sinw/Xi cosa;^a:;2+ 
cos ujfXi sinwja;2) 
sinw/xs (cosw/Xi cosa;/a;2 — 
sinw/Xi sinLLi/a;2) . 

The other two bispectra are not of the same form, since 
their arguments are non-hnear functions, i.e. Inxt for 
resonant non-Gaussianities and (xj+i + Xj+2)/xj for for 
non-BD modifications and these can be made of the form 
above by expanding, using the constructed Fourier modes 
J-fi ■ Given the form of the first argument you expect only 
a limited number of modes to significantly contribute, for 
example those modes that have equal mode number in 
the directions xi, X2 and X3 (you should think about this 
expansion as a series around the point Xt , see table |ll| . 
For the second argument you expect more modes to mat- 
ter, since the arguments depend on all three directions in- 
dependently. Consequently the weights Q!„ are expected 
to be close to zero for many n when expanding resonant 
non-Gaussianities, while for the non-BD scenario they 
should all matter to some extend (and obviously more 
modes will be important for large uj^). 



FIG. 6: The one dimensional orthonormal Fourier functions 
fn{xi) within the tetrahedral domain for the first 11 modes. 

refer to the three dimensional orthonormalized modes as 
J-„ and the corresponding weights as Q!„. 

N 

S{xi,X2,X3) ~ ^TZe{anJ^ri{xi,X2,X3)) . 
n=0 

If S would have been complex, one should add 
iXm (a„J-"„) in order to take this into account. The coef- 
ficients Un can be computed by taking the inner product 
(eq. ( [T5| ) between the original shape function (bispec- 
trum) and the various mode functions J-m i.e. 

a„ = / dxidx2dx3S{xi,X2,X3)J^*{xi,X2,X3). 



1. Feature Bispectrum 

For the feature bispectrum we do not necessarily have 
to consider the Fourier expansion^ since that spectrum 



We will later show that we would also find it when wc would 
search for resonant bispectra in the data, as the weights peak at 
almost the same mode numbers. 



2. Resonant Bispectrum 

We have computed the correlation for the axion bispec- 
trum with the Fourier expansion for frequency ranges of 
= 20 — 80 up to 82 modes (figure [?]). As expected, we 
see that there are only a few modes that give significant 



1.0 




FIG. 7: Correlation between sin(tJr lna;t) for frequencies ujr = 
20, 40, 60 and 80. Compared with the polynomial mode ex- 
pansion we reach similar correlation using about 5 times less 
modes. Also note that the increase of correlation is some- 
what discreet, indicating that we might need only a fraction 
of these modes to reconstruct the original spectrum. We will 
discuss this observation in the next section. 
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FIG. 8: The correlation between the DBI bispectrum and 
both polynomial and Fourier expansion as a function of the 
maximum number of modes. As expected, both expansions 
start out equally (the zero mode of the real part of the Fourier 
expansion is equivalent to the zero mode of the power law 
expansion), but while the power law reaches a correlation of 
> 99% after just 5 modes, the Fourier remains stuck at 97%. 



contribution to the correlation, while most modes give 
only very little contribution and are not important for 
the expansion. We will discuss this fact in the context 
of CMB data mode extraction in the section 4. Given 
that the allowed range of frequencies 20 < ^ 10^ this 
expansion is actually reasonable for the lower frequencies 
and the number of modes necessary to establish similar 
correlation as the polynomial expansion is reduced by a 
factor 5. 



3. Non-BD Bispectrum 

As for the polynomial basis expansion, the presence 
of a large number of features in the non-BD bispectrum 
does not allow for a fast reconstruction of the spectrum. 
In fact, expansion in the Fourier basis requires even more 
modes compared to the polynomial basis, reaching only 
~ 20% correlation after 82 modes with — 20. We 
also find that Uf — actually reaches a slightly larger 
correlation, although this seems mostly due to a rela- 
tively large correlation with the zero order (n = 0) mode. 
Most likely this is caused by the fastest oscillating part of 
the spectrum which, in combination with numerics, could 
add constant power. We did observe something similar 
in figure [5] for polynomial modes where the zero mode 
causes the correlation of the non-BD bispectrum recon- 
struction with uj^ — 40 to be better initially compared to 
bispectrum expansion with ujf = 20. 

In most realistic scenarios uj^ > 100 (otherwise your 
effective field theory approach breaks down) and there- 
fore both polynomial expansion and Fourier expansion 
fail to reconstruct this bispectrum effectively. The pos- 
sible explanation why Fourier expansion is even worse 



than polynomial expansion for this type of bispectrum, 
seems to be related to the rapid change in frequency in a 
fixed direction. Fourier expansion is optimized for scale 
invariant frequencies. The polynomial expansion is sim- 
ply optimized in reproducing as many different shapes 
as possible, explaining the observation that it is able to 
slowly increase correlation with the addition of modes 
while Fourier expansion seems to converge. Given the 
large enhancement of the amplitude /^'f^ (which scales 
as w^), one might still be able to extract some informa- 
tion from that data even with such small correlations. 

Another possibility is that once non-BD bispectrum 
is projected onto multipole space one might establish 
a larger correlation with fewer (multipole) modes. The 
projection has the tendency to wash out small features 
(hence the weight of l/Zcj in the correlator.). We hope to 
report on this in the future. 



4- Toy Spectra 

To investigate the power of the Fourier expansion for 
oscillatory bispectra we have also tried to fit three toy- 
model shapes moving in different direction through co- 
moving momentum space 
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(19) 



We find again that for such a shapes the correlation 
increases about 5 times faster compared to polynomial 
mode expansion with the same frequency. The corre- 
lation as a function of mode numbers for Fi and F2 are 
shown figure |9] We will discuss the weights of these mod- 
els in the next section. Note that Fi is already of the fac- 
torized form, however here we simply aim at showing the 
effectiveness of Fourier expansion. We want to emphasize 
that these spectra are not based on any physical model, 
but simply show that in general oscillatory spectra are 
better fitted using a Fourier basis. 



5. Smooth Spectra 

Although the Fourier expansion seems to work well 
for resonant non-Gaussianities and the toy-spectra, com- 
pared to polynomial expansion we confirm that Fourier 
expansion is not as effective: it is easier to gain fast con- 
vergence with a limited number of modes for most os- 
cillating bispectra, but it is difficult to get correlation 
beyond 0.97 for smooth bispectra. This is probably due 
to overshooting at the boundaries as discussed in [53]. 
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FIG. 9: The correlation as a function of mode number for two out of three toy-spectra Fi(left) and _F2 (right) in eq. (191. In 
both cases Fourier expansion (dashed) leads to faster convergence compared to polynomial expansion (solid). 



We explicitly show this in figure |8] where we compare ex- 
pansion of the 'smooth' DBl inflation bispectrum (which 
is very similar to equilateral), using Fourier modes and 
polynomial modes. 

We conclude that Fourier expansion is a viable alter- 
native for polynomial expansion in the case of oscilla- 
tory bispectra with relatively large frequencies. Using the 
Fourier expansion we can achieve factorizabilty of various 
oscillating bispectra with significantly less modes com- 
pared to polynomial expansion. For frequencies w 3> 50 
polynomial and Fourier expansion are both unable to 
reconstruct the original spectrum with a small number 
of modes. In order to reconstruct models with such 
large frequencies, one should look for alternative meth- 
ods. However, constraining these models with only lim- 
ited number of modes seems to be a practical possibility. 
This will be topic of the next section. 



out of three directions have one and two maximal orders 
less than the third, i.e. in mode number n = 23 two 
directions are maximally quadratic and the third is max- 
imally cubic. As we already argued the location of these 
peaks makes sense, since the resonant model is a func- 
tion of kt (or Xt), which is the sum of the three comoving 
momenta. Effectively this shape is orientated in the kt 
direction. One could only try to expand the spectrum 
only in those modes, which could significantly reduce the 
number of modes necessary. Since the important modes 
seem to be related to the direction of propagation of the 
oscillation, we find that this conclusion is independent of 
the phase. In other words, only the value of the weights 
will differ, not the mode numbers that are relevant for 
the expansion. 

Secondly, from an observational point of view, given 



IV. DISCUSSION 

Even though the expansion of the oscillatory primor- 
dial bispectra becomes unavailing for really large frequen- 
cies, there are a number of interesting observations which 
could make constraining and expanding oscillating bis- 
pectra much more viable than presently argued. First 
of all, as predicted, the expansion in mode functions of 
the resonant bispectrum has a very discrete character; 
basically if you consider fig. [7] only few modes actually 
contribute significantly to the convergence of the corre- 
lation. In fig. 11 we show the various weights (|5„|) 
as a function of mode number (as well as for Fi and 
F3 (not shown)). We can trace back the corresponding 
mode numbers in table |Tll For instance there is a clear 
peak at n = 16, which correspond to all directions being 
maximally of quadratic order, and n = 41 with all direc- 
tions being maximally of cubic order. Other peaks (e.g. 
n = 21, 32 and 53) correspond to the modes in which two 
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FIG. 10: The weights derived for 65 modes for both cosujr In kt 
(dashed) with cj,. = 50 and sin ujfkt (solid) with oj/ = 20 
showing that these both peak for similar mode numbers. Al- 
though distinguishing between these would be quite hard, it 
seems that for the feature bispectrum the values of the weights 
fin are peaked sharper. 
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FIG. 11: The weights ^QnQ* for a resonant bispectrum as a function of mode number for various frequencies. It is clear 
that only limited number of modes are valuable in the reconstruction of the original spectrum via mode expansion. From an 
observational point of view this is very convenient as it would require the measurement of only a limited number of Fourier 
modes to learn about oscillations in the primordial bispectrum. 



the discreteness of the correlation it is (obviously) not 
necessary to constrain all mode functions in the CMB 
data to get an indication of there is an oscillatory three 
point signal and what the possible frequency of this signal 
might be. For resonant non-Gaussianities we only need 
to consider those modes that have a significant weight 
a, and the measured value of the weights would be a di- 
rect measure of the frequency. If one could extract the 
multipole projected Fourier modes that are responsible 
for most of the weight, this could in principle provide 
signatures of primordial bispectra with frequencies much 
larger than cj^ = 80. Measuring modes up to e.g. n — 100 
would not only provide information about the frequency 
of the signal, but could also hint on the type of primordial 
bispectrum. The distinction between the feature bispec- 
trum Q and the resonant bispectrum ([9| would be more 
difficult, since the values of the weights peak at similar 
mode numbers although we have found that expanding 
the feature bispectrum in the constructed Fourier basis 
(instead of the simple trigonometric expansion discussed 
in section [ill B 1 ) could still be used to discriminate be- 
tween the two signals (see figure 10). 

To emphasize the ability to extract information on the 
primordial shape solely from the modes that are impor- 
tant, we have investigated three toy-model shapes of eq. 



We have computed the Fourier weights for two 
different frequencies in figure 12 As expected, Fi has 
weights that peak when only one comoving momentum in 
in k space is non-zero is, i.e. it peaks at the modes where 
one momentum oscillates and the other two momenta are 
constant (see|ll]). The obvious reason is that each term in 
Fi depends on one comoving momentum variable only, 
implying that there should be no cross terms in the ex- 
pansion. For F2 we find that many more modes are rele- 
vant, which makes perfect sense given that the argument 
in the sine depends on all three vectors in comoving mo- 
mentum space. For F^ however the argument effectively 
only depends on two comoving momenta, therefore the 
relevant mode functions (the ones with the largest |q;|) 
are the ones that have similar frequency in two momen- 
tum vectors and are constant in the third. 

In this paper we have only discussed mode functions in 
momentum space, and one either has to construct similar 
Fourier modes in multipole space or project these modes 
forward using the transfer function |23j , and use these to 
expand a late-time oscillatory bispectrum, and see if we 
get similar results in terms of mode number sensitivity. 
One expects that after projection the transfer function 
has caused some smoothing of the signal, which could 
render a Fourier basis less effective. On the other hand. 
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FIG. 12: The weights -^Q„aJ for the three toy-spectra as a function of mode number for various frequencies. For the top 
spectrum the most relevant modes are those that have a maximal frequency of ~ 27m for one comoving vector and are constant 
for the other two. For the second toy-spectrum (middle) the symmetry of the argument makes all weights relevant (decreasing 
as a function of mode number). The third example (bottom) has modes dominating of which two comoving momentum 
vectors have non-zero (equal) frequency and the third one is constant. For example, the lower-right bottom example has jasTi 
dominating, which corresponds to {p, r, s} = {0, 5, 5}. 



intuitively it seems perfectly reasonable that a Fourier 
basis should be much more efficient in reconstructing os- 
cillatory bispectra from the data. In addition, the effects 
of the transfer function on the correlation in I space can 
be examined by choosing the 1/kt weight Wk in the pri- 
mordial correlation function. We have found that our 
results were only marginally affected when including this 
weight factor and therefore we expect that Fourier mode 
expansion should be equally efficient in multipole space. 
To make sure that this is actually true, we should com- 
pute the projection of several oscillatory bispectra and 



construct a orthonormal Fourier basis in multipole space. 
We hope to report on this in the near future. 



V. CONCLUSIONS 

We have investigated the viability of mode expansion 
for bispectra that contain oscillations. The motivation 
for investigating such features and their mode expansion, 
is that recently it has been shown that several scenarios 
or mechanisms can produce such features not only in the 
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power spectrum, but also in the bispectrum. The appear- 
ance of oscillations in the bispectrum makes comparison 
with existing bispectral constraints, based on smooth bis- 
pectra, very inefficient and there exists substantial room 
for improvement. In order to constrain oscillatory bispec- 
tra from the data, a logical first step is to factorize the 
bispectrum in order to efficiently compute its multipole 
counterpart. Polynomial expansion has been proposed 
to achieve factorization of a given theoretical bispectrum 
and we have investigated this for three different models. 
As expected, the larger the frequency of the primordial 
bispectrum, the more modes it requires to establish a 
reasonable approximation of the original spectrum. In 
the case of a feature in the primordial potential polyno- 
mial mode expansion might still be useful, at least for 
features at high multipoles (resulting in rather small fre- 
quencies in comoving momentum space). In fact, during 
the finalization of this paper the authors of |31j have 
considered a feature bispectrum and extracted 31 poly- 
nomial modes in the data, which allowed them to inves- 
tigate late time bispectra with a maximal frequency of 
ujf = 5—10 (in comoving momentum space). They did 
not find 3(7 evidence for non-zero non-Gaussianity. The 
other two example bispectra typically have a lot more 
oscillations within the tetrahedral domain, resulting in 
many modes necessary to realize an acceptable correla- 
tion. Fortunately, both the resonant and non-BD bispec- 
trum have an amplitude that scales with the frequency. 
Therefore, a small improvement in correlation could lead 
to a significant improvement in the ability to constrain 
the model by measuring these modes in the data and 
reconstructing the primordial signal. 

Complementarily, we have proposed a different basis 
expansion, based on Fourier functions instead of poly- 
nomials. This still leads to the necessary computational 
reduction one is after and therefore is a perfectly valid 
alternative. Such expansion is more relevant for resonant 
and non-BD scenario, since the feature bispectrum can 
already be transformed into Fourier modes analytically, 
using identities. We have shown that Fourier modes are 
much more efficient for the resonant bispectrum, reduc- 
ing the number of modes necessary to establish the same 
correlation as polynomial modes by at least a factor of 
5. For the non-BD bispectrum both Fourier expansion 
and polynomial expansion are difficult. Correlation in- 
creases fast with the addition of modes, but quickly con- 
verges to a fixed value, where the fixed value decreases 
a function of frequency. We believe that this is due to 
the exact form of the bispectrum, which has many small 
features near the edges of the tetrahedral domain. One 
might hope that some of these very small features are 
washed out when you compute the multipole equivalent. 



although that would be very time consuming since the 
non-BD shape is not of the factorized form. We hope to 
investigate this in a future attempt. In addition we have 
investigated three toy-spectra, not based on any partic- 
ular model, which have a different oscillating orientation 
compared to the three theoretical models. Expanding 
these in Fourier modes show similar improvement com- 
pared to polynomial expansion as the resonant bispec- 
trum. In general, we therefore belief that Fourier expan- 
sion is much more effective in the expansion of oscillatory 
spectra compared to polynomial basis expansion. 

From an observational stand point, it seems that for 
resonant inflation only a limited number of modes con- 
tribute significantly in reproducing the original bispec- 
trum. This allows us to consider only those modes that 
contribute substantially. This holds independent of the 
phase and frequency of the signal and is due to the spe- 
cific form of this bispectrum, which oscillates (primarily) 
in the kt direction. Because the modes that are impor- 
tant for the reconstruction of the original bispectrum are 
independent of the frequency, this also implies that when 
one would observe these modes in the data one could in 
fact flnd evidence for much larger frequencies than dis- 
cussed here, simply because for larger frequencies these 
modes will also matter but their respective weight will 
be smaller. Despite the fact that we could not optimally 
expand the non-BD bispectrum using Fourier modes, we 
did look into the three toy-sepctra. We found that other 
modes are important. Moreover, the modes that are im- 
portant directly represent the orientation of the oscillat- 
ing spectrum and could therefore discriminate between 
different bispectra quite effectively. If this conclusion 
holds after forward projection into multipole space, mea- 
suring a number of Fourier mode functions in the CMB 
data would present an efficient way of deducing whether 
oscillations are present in the data and could give both 
an indication of the frequency and the shape of the pri- 
mordial bispectrum. 
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